The Vogel-Fulcher-Tammann law in the elastic theory of glass transition 
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We propose that the origin of the Vogel-Fulcher-Tammann law is the increase of the range of 
elastic interaction between local relaxation events in a liquid. In this picture, we discuss the origin 
of cooperativity of relaxation, the absence of divergence of relaxation time at a finite temperature 
and the crossover to a more Arrhenius behaviour at low temperature. 
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The transition of a liquid into a glass on lowering the 
temperature may appear conceptually simple, yet this 
phenomenon has turned out to be one of the most difficult 
and controversial problems in condensed matter physics, 
the problem of the glass transition [H, Q . At high tem- 
perature, relaxation time r of a liquid follows Arrhenius 
dependence. On lowering the temperature, r almost uni- 
versally deviates from Arrhenius dependence, and follows 
the Vogel-Fulcher-Tammann (VFT) law: 



r exp 



.4 



T-T 
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where A and To are constants. The origin of the VFT 
law is the main open question in the field of the glass 
transition [l], 0] • 

A related open question follows from the form of the 
VFT law, namely what happens at To. Because r for- 
mally diverges at To, several models have suggested that 
a phase transition from a liquid to a glass phase can ex- 
ist P, 0] . Because the divergence is not observed in an 
experiment, it was proposed that the phase transition is 
avoided due to sluggish dynamics when r exceeds exper- 
imental time scale. However, the nature of the phase 
transition and the second phase is not clear, which con- 
tinues to fuel the current debate [H, Q ■ Interestingly, the 
VFT law changes to a more Arrhenius form at low tem- 
perature, pushing the divergence temperature down 
The origin of this crossover is not understood. 

Another related problem is the physical origin of "co- 
operativity" . The notion of cooperativity of molecular 
motion, which sets in a liquid as temperature is lowered, 
was introduced and intensely discussed in several popular 
theories of the glass transition. These theories are based 
on the assumption that "cooperatively rearranging re- 
gions" , "domains" or "clusters" exist in a liquid, in which 
atoms move in some concerted way that distinguishes 
these regions from their surroundings [H, [3, 0, a [E 13 • 
The physical origin of cooperativity is not understood, 
nor is the nature of concerted motion. 

A glass is different from a liquid by virtue of its ability 
to support shear stress. This suggests that the change 
of stress relaxation mechanism in a liquid on lowering 
the temperature is central to the glass transition pro- 
cess, yet stress relaxation is not discussed in popular 
glass transition theories, including entropy, free-volume, 




FIG. 1: An example of a concordant local relaxation event. 
Solid and dashed lines around the shaded atoms correspond 
to initial and final positions of a rearrangement, respectively. 
Arrows show the direction of external stress. 



energy landscape and other approaches Q. 

In this paper, we discuss how stress relaxation in a 
liquid changes with temperature. We propose that the 
origin of the VFT law is the increase of the range of elas- 
tic interaction between local relaxation events. In this 
theory, we also discuss the origin of cooperativity of re- 
laxation, the absence of divergence of r at a finite temper- 
ature and the crossover to a more Arrhenius behaviour 
at low temperature. 

Relaxation and flow in a liquid proceed by elemen- 
tary local structural rearrangements, during which atoms 
jump out of their cages. We call these rearrangements lo- 
cal relaxation events (LREs). Because the divergence of 
the elastic field due to a LRE is zero, a LRE is not ac- 
companied by compression of the surrounding liquid, and 
can be viewed, in a simple model, as a pure shear event 
[2J. Therefore, in discussing how LREs interact elasti- 
cally, we consider shear LREs. A typical shear relaxation 
event is shown in Figure 1 (term "concordant" in the fig- 
ure caption is not important here, and will be explained 
later). The accompanied structural rearrangement pro- 
duces elastic shear stress which propagates through the 
system and affects the relaxation of other events. The 
important question here is how does this stress affect re- 
laxation of other LREs in the liquid? 

Lets consider how the changes of stresses due to re- 
mote shear LREs affect a given local relaxing region, 
shown in the centre in Figure 2. Relaxation of the cen- 
tral event involves deformation of the "cage" around the 
jumping atom (see Figure 1), and therefore depends on 
the stresses that propagate from the remote LREs to the 
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teratomic separation of about 1 A and To the oscillation 
period, or inverse of Debye frequency (tq ~ 0.1 ps), 



FIG. 2: Illustration of the elastic interaction between local re- 
laxation events. This interaction takes place within the range 
d e i from the central relaxing regions. Shaded and open circles 
represent local relaxing regions inside and outside, respec- 
tively, of the interaction sphere. 



centre. A remote shear LRE, similar to the one shown in 
Figure 1, creates elastic shear waves, which include waves 
of high frequency This is because the deformation, asso- 
ciated with a LRE, creates a wave with a length compa- 
rable to interatomic separations (see Figure 1), and hence 
with a frequency on the order of the Debye frequency. At 
high frequency u> > 1/r, a liquid supports propagating 
shear waves [8[ , which propagate stress and its variations 
from remote LREs to the central point. If r is macro- 
scopically defined as the time of decay of shear stress in 
a liquid [1, Q , d c \ = ct gives the length of this decay, 
where c is the speed of sound. Here, d c i gives an estima- 
tion of the maximal range over which shear stress decays 
in a liquid. At the microscopic level, the relevance of 
d c \ = ct is as follows. A high-frequency shear wave origi- 
nating from a LRE propagates stress until a remote LRE 
takes place at the front of the wave, at which point the 
wave front is absorbed by the remote LRE. Suppose this 
happens at distance d e \ from the original LRE. d c \ can 
be calculated from the condition of equality of the wave 
travel time, d e i/c, and the time at which the remote LRE 
takes place at point d c \. The latter time is given by r, 
because microscopically, r is defined as the average time 
between two consecutive LREs at one point in space Q, 
and we obtain d c \ — ct as before. 

Therefore, d c \ defines the maximal distance over which 
the central LRE is affected by elastic shear stresses due 
to other LREs in a liquid (see Figure 2). For this reason, 
d c \ can be called the liquid elasticity length. Note that 
relaxation of the central event is affected by all those 
stresses that have enough time to propagate to the centre. 
Because it takes time r for the central event to relax, 
its relaxation is affected by the stresses from all LREs 
located distance ct away. After time r, the central event 
relaxes, and the process repeats. Therefore, the definition 
d c \ — ct is self-consistent. 

Because c is on the order of q/tq, where a is the in- 
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TO 



(2) 



On lowering the temperature, r increases as t = 
Tnex-p(V/kT), where V is the activation barrier of a LRE 
[8( (here, V can be temperature-dependent). According 
to Eq. (2), this increases d c \ and the number of LREs 
that elastically interact with a given event. We propose 
that this is the key to the super- Arrhenius relaxation. 

Before discussing the VFT law itself, we note that Eq. 
(2) immediately gives the crossover from non-cooperative 
to cooperative relaxation. When, at high temperature, 
t w To, d e \ ~ a (see Eq. (2)), and d c \ < d m , where 
d m is the distance between neighbouring LREs of about 
10 A {d m is the distance between the centres of neigh- 
bouring molecular cages). This means that LREs do not 
elastically interact. As r increases on lowering the tem- 
perature, d i > d m becomes true. At this point, LREs 
are no longer independent, because relaxation of a LRE 
is affected by elastic stresses from other events. This 
discussion, therefore, clarifies the physical origin of co- 
operativity. Here, we do not need to assume or postu- 
late cooperativity of relaxation as in the previous work 
[E 0, A H, B 0- I n this picture, relaxation is "coop- 
erative" in the general sense that LREs are not inde- 
pendent, but the origin of this cooperativity is the usual 
elastic interaction. We have recently shown how this in- 
teraction gives rise to stretched-exponential relaxation 
(SER), a universal feature of supercooled liquids [Tpj ] . 
The crossover from exponential relaxation to SER takes 
place when d e \ — d m . According to Eq. (2), r at the 
crossover, t c , is a universal value: t c — Tod m 

/a. This 

gives t c of about 1 ps, consistent with the numerous ex- 
periments (ill [i2T|. 

In order to derive the VFT law, we recall the previous 
discussion that V is given by the elastic shear energy of a 
liquid around a LRE [3, [l^, 14 1 . The energy needed for an 
atom to escape its cage at the constant volume is very 
large because of the strong short-range interatomic re- 
pulsions, hence it is more energetically favourable for the 
cage to expand, reducing the energy needed for escape. 
Such an expansion elastically deforms the surrounding 
liquid, hence V is given by the work of the elastic force 
needed to deform the liquid around a LRE. Because this 
deformation does not result in the compression of the sur- 
rounding liquid (for the displacement field u created by 
an expanding sphere, div(u) = 0), V is given by the back- 
ground shear energy of the liquid. This was confirmed by 
the experiments showing that V increases with the liquid 
shear energy 
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We now recall the previous discussion of how LREs re- 
distribute external stress. In discussing creep, Orowan 
introduced "condordant" LREs [15| . A concordant shear 
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LRE is accompanied by a strain in the direction agreeing 
with the applied external stress, and reduces the local 
stress and energy (see Figure 1). In order to counter- 
balance this decrease, other local regions in a system sup- 
port more stress [lEj]. Goldstein applied the same argu- 
ment to a viscous liquid under external stress Con- 
sider that this stress is counterbalanced by stresses sup- 
ported by local regions. Because a local region supports 
less stress after a concordant LRE than before, other lo- 
cal regions in the liquid should support more stress after 
that event than before in order to counter-balance jl6| . 

Lets consider a liquid perturbed by a pulse of an ex- 
ternal field. At time zero, shear stresses supported by 
local regions counterbalance external shear stress. As re- 
laxation proceeds, each concordant shear LRE reduces 
stress locally, until the external stress is relaxed by a cer- 
tain number of LREs N. When this process is complete, 
the liquid relaxes to equilibrium. At times smaller than 
L/c, where L is the system size, the external stress can 
be considered constant, and the stress redistribution ar- 
gument of Orowan-Goldstein applies. Alternatively, we 
can consider an external stress constantly compensating 
for the decreases of local stresses. In the resulting steady 
flow, r is the time needed to relax an increment of exter- 
nal perturbation, and can be viewed as the time of the 
liquid's retardation behind the external field. Let n be 
the current number of LREs, such that n — > N. If Ap is 
the increase of shear stress on the liquid around a current 
local region that arises from the need to counter-balance 
the decreases of stresses due to previous remote concor- 
dant LREs, Ap increases with n. The increase of Ap 
consistently increases the elastic strain in the direction 
of external shear stress, increasing the background shear 
energy of the liquid around the current local region. As 
discussed above, V for the current LRE increases as a 
result. The increase of V, AV, due to Ap is the work 
J Apdq. If q a is the activation volume [14|, AV — Apq a , 
and V = Vo + q a Ap, where Vq is the high-temperature 
activation barrier. Because Ap increases with n, V also 
increases with n. This gives the elastic feed-forward in- 
teraction mechanism for LREs, which sets SER (Toj . 

To calculate V as a function of d e i, lets consider the last 
LRE that relaxes an increment of external shear stress to 
be in the centre of a sphere of radius d e \ (see Figure 2). 
As relaxation proceeds, the shear stress on the central re- 
gion increases in order to counterbalance stress decreases 
due to previous remote concordant LREs. Importantly, 
because this mechanism operates in the range set by d e \ 
and because d e i increases on lowering the temperature 
(see Eq. (2)), stresses from an increasing number of re- 
mote LREs need to be counterbalanced by the central re- 
gion. It is also important to note that all stresses within 
a distance d e \ — ct have enough time to propagate to the 
centre and affect relaxation of the central event (recall 
self-consistency in definition of d e \). 

Let Apj(O) be the reduction of local stress due to a re- 



mote concordant LRE i. Api decays with distance, hence 
we denote Apj(r) as its value at the centre in Figure 2. 
The increase of stress on the central rearranging region, 
Ap, can be calculated as 



Ap = p J Aitr 2 Api(r)dr 

d /2 



(3) 



where p is the density of local rearranging regions and 
do is on the order of the size of a relaxing region (in 
Figure 1, do > 3a). Note that in Eq. (3), d e l is the 
upper integration limit. In what follows, we assume, for 
simplicity, that Api(0) are constant, Api(O) = Apo- 

In an elastic medium, stresses decay as Ap(r) oc 
1/r 3 @. Because Ap(r) = Ap at do/2, Ap(r) = 
Apo(do/2r) 3 . Integration of Eq. (3), together with 
V = Vo + q a Ap from the discussion above, gives 

V = V + Tr/2pq a Ap dl ln{2d el /d ) (4) 
Using t — tq exp(V/kT) in Eq. (2), we obtain 
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Eqs. (4) and (5) define V in a self- consistent way. 
Eliminating d e i from the two equations, we find: 



V = 



AT 

T - To 



(6) 



where A = Vo + Tr/2pq a Apod i \n(2a/do) and kT 
Tr/2pq a Ap d i . 

From Eq. (6), the VFT law follows: 



T exp 
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In this picture, the super-Arrhcnius behaviour is re- 
lated to the increase of d e \ (see Eq. (4)). The transition 
from the VFT law to the Arrhenius form of r takes place 
in the limit of small d c \ at high temperature. In this 
case, the upper and lower integration limits in Eq. (3) 
coincide, giving Ap = 0, V = Vo and r = r exp(Vb/fcT). 

In the proposed theory of the glass transition, the on- 
going controversy [l|, [H, [l8| regarding the divergence and 
possible phase transition at To is readily reconciled. The 
divergence at To can not exist for the following reason. 
From Eqs. (5,6), we find 



dei = a exp 



T-Tn 



(8) 
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When T approaches To, d e \ diverges, and quickly exceeds 
any finite size of the system L. When d e \ > L, all LREs 
in the system elastically interact, and there is no room 
for the increase of V by way of increasing d c \. The upper 
limit of integral (3) becomes d e \ = L, giving temperature- 
independent V oc ln(L) (see Eq. (4)). Further decrease 
of temperature has a weaker effect on V, and can be due 
to, e.g., density increase, but not to the increase of d c \ 
(the density-related contribution to V does not depend 
on d e } ° r L). As a result, the behaviour of r tends to 
Arrhenius, pushing the divergence to zero temperature. 

d e \ exceeds the experimental value of L above T g : if 
T (T g ) = 10 3 sec, d ci {T g ) = 10 3 km, according to Eq. (2). 
Hence our theory predicts the crossover from the VFT 
law to a more Arrhenius behaviour at low temperature, 
as is seen in the experiments 0. According to Eq. (2), 
r at the crossover is r = tqL/ a. If a typical value of L is 
1 mm, r at the crossover is 10 -6 sec, consistent with the 
experimental results [l7j ]. 

We note here that d e i vastly exceeds the size of "co- 
operatively rearranging regions" (CRR), which is several 
nm at T g (for review, see, e.g^ Ref. The physical 

picture of CRR is not clear [2j. It is possible that the 
observed nm scale of CRR is set by the distance beyond 
which the elastic strains from LREs decay to the values 
undistinguishable from thermal fluctuations. 

d c \ gives an insight into the origin of liquid fragility 
[l8j |. According to Eq. (4), as long as at high temper- 
ature d c \ < L, lowering the temperature increases V, 
resulting in a fragile behaviour. If, on the other hand, 
del > £ at high temperature already, further decrease of 
temperature has a weaker effect on V , giving weak super- 
Arrhenius behaviour. Experimentally, for many systems 
the studied range of temperatures varies from about 2T g 
and T g hence we consider the increase of d e i from 
high temperature Th — 2T g to T g . Take, for example, two 
systems on the stronger side of fragility plots, BeF2 and 
SiC>2. From the experimental values of Vh/kT g (Vh is the 
activation barrier at the highest measured temperature), 
we find V h JkT h = 24 and 19.6 for BeF 2 and Si0 2 , re- 
spectively [19(. According to Eq. (5), this gives d e i = 2.6 
m and 33 mm at Th for the two systems. Because a typ- 
ical experimental value of L is on the order of 1 mm, our 
theory correctly predicts that these systems should be 
on the strong end of fragility plots. For two fragile sys- 
tems, toluene and propylene carbonate, Vh/kTh = 3.34 
and 5.75, giving d e i — 28 and 314 A at Th, respectively. 
This is much smaller than L, hence our theory predicts 
that these systems should be fragile, as is seen experi- 
mentally. An interesting prediction from this picture is 
that strong systems will show increased fragility at high 
temperature when d c i < L (note that strong systems have 
been measured at relatively low temperature only [181]). 

Before concluding, we note that we discussed a typical 
experimental setup, in which a liquid is perturbed and 
r is measured as the time of return to equilibrium. All 



above results remain the same in the equilibrium case 
as well, when thermally fluctuating LREs interact via 
induced elastic stresses in the range set by d e i (iol |. 

In summary, we proposed that the origin of the Vogel- 
Fulcher-Tammann law is the increase of the range of elas- 
tic interaction between local relaxation events in a liquid. 
In this picture, we discussed the origin of cooperativity of 
relaxation, the absence of divergence of relaxation time 
at a finite temperature and the crossover to an Arrhenius 
behaviour at low temperature. 

We suggest that the proposed theory is applicable to 
other systems, in which local rearrangements interact 
via the fields they induce. This includes a wide range 
of phenomena, for example, relaxation in spin glasses. 
Here, the same universal relaxation effects as in struc- 
tural glasses are observed, including the VFT law, coop- 
erativity, SER and other phenomena. 
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